A simple nonlinear equation for structural relaxation in glasses 
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A wide range of glassy and disordered materials exhibit complex, non-exponential, structural 
relaxation (aging). We propose a simple nonlinear rate equation 8 = a [1— exp (bS)], where 8 is 
the normalized deviation of a macroscopic variable from its equilibrium value, to describe glassy 
relaxation. Analysis of extensive experimental data shows that this equation quantitatively captures 
structural relaxation, where a and b are both temperature-, and more importantly, history-dependent 
parameters. This analysis explicitly demonstrates that structural relaxation cannot be accurately 
described by a single non-equilibrium variable. Relaxation rates extracted from the data imply the 
existence of cooperative rearrangements on a super- molecular scale. 



Glass-forming materials exhibit rapid increase in relax- 
ation timescales when going through their glass temper- 
ature T g When external conditions change abruptly, 
the observation of the full relaxation to equilibrium be- 
comes exceedingly difficult, except for a narrow range of 
temperatures near T g . The relaxation is characteristi- 
cally non-exponential and spans several orders of magni- 
tude in time. In many cases, the relaxation is logarith- 
mic. This behavior is paralleled in a broad range of dis- 
ordered systems: compaction of granular materials 0,0], 
crumpling of thin sheets [H, Q , aging of contact area in 
dry friction [(|, aging of conductivity in electron glasses 
7, |8| and mechanical relaxation of star polymers in gels 
II llfj . This apparently wide-spread behavior might sug- 
gest a generic origin of slow glassy relaxations (aging). 

Glassy relaxation is typically probed by tracking the 
time evolution of a macroscopic quantity, e.g. the volume 
or the enthalpy of a sample, in response to an abrupt 
change in an externally controlled variable, e.g. the tem- 
perature. In the latter case, when an initial temperature 
To is rapidly changed to T, in the vicinity of the glass 
temperature T g , various degrees of freedom of the glass 
respond differently. The vibrational degrees of freedom 
quickly equilibrate at the new temperature T. The struc- 
tural degrees of freedom, however, carry long-time "mem- 
ory" of the original state at To and fall out-of-equilibrium 
with T. It is the out-of-equilibrium dynamics of the 
structural degrees of freedom towards a new equilibrium 
at T that is at the heart of "structural relaxation" . 

Structural relaxation is conventionally interpreted in 
terms of the Tool-Narayanaswamy-Moynihan (TNM) 
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or its equivalent KAHR [12j, phenomenological 4- 
parameter models. The main assumption in these models 
is that during relaxation, a single dynamical variable is 
sufficient for describing the non-equilibrium state of the 
glass. In the case of the TNM model this variable is a 
"fictive temperature" , defined to be a linear function of 
the probed property, e.g. volume. 

However, it is known that these models do not de- 
scribe experimental data accurately [l| , and actually fail 



relaxation (i.e. the Kovacs memory effect [1J]) cannot 
be described by a single non-equilibrium state variable 
15 , 1(| , some recent works have suggested that this might 



be possible for monqtonic relaxations of the type consid- 
ered here 



15, 17, 1 



In spite of the seemingly universal nature of structural 
relaxation in glasses, as well as its great scientific and 
technological importance, a theoretical understanding of 
it is still missing. In this Letter, we propose a simple, 
analytically solvable, nonlinear rate equation for describ- 
ing structural relaxation (aging) in glasses. The pro- 
posed equation is shown to quantitatively capture exten- 
sive experimental measurements on volume relaxation, 
yet it explicitly demonstrates the inadequacy of a single 
non-equilibrium description of glassy relaxation. This 
analysis allows us to shed light on some basic proper- 
ties of structural relaxation in these systems, including 
estimates of activation energy barriers and volumes. 

Volume relaxation of glassy materials is usually stud- 
ied using mercury dilatometry fl9j | . The classical exper- 
iments in volume relaxation were performed by Kovacs 
more than five decades ago [l4T |. These measurements, 
which in extreme cases reached months, are routinely 
used and provide a standard testing ground for models. 
In these so-called down-jump experiments, a glass sam- 
ple is rapidly quenched from equilibrium at To to a lower 
temperature T. Measurements begin at ij, the character- 
istic time it takes for the vibrational degrees of freedom 
to thermalize. 

The first question we raise is whether structural re- 
laxation can be properly described by a single non- 
equilibrium variable model. To address this question, 
we denote such a variable by 5(t) and note that a sin- 
gle variable description means that the rate of relax- 
ation is uniquely determined by the instantaneous value 
of the variable, i.e. 5(t) = r[6(t)], where r[-] is some func- 
tional. Therefore, in the framework of such models there 



exists a function g(t) such that S(t) = g(t 



to account for thermal history dependence j 1 31 ] - In ad 



dition, while it has been recognized that non-monotonic 



where Sq = 5(0) is the initial condition, and hence a 
set of measurements which differ only in initial condi- 
tion could be time-shifted in such a way that all curves 
would collapse on a single master curve. An example of 
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such a set is shown in Fig. QJa) with measurements dig- 
itized from Kovacs' original work 14 1. given in terms of 
8= (V — Voo ) /Voo , the normalized deviation of the volume 
from its asymptotically stable value V^. Here the tem- 
perature T = 30°C is the same for all measurements, and 
only the initial state of the system is changed by quench- 
ing from different initial temperatures To- In Fig- Hlb)' 
we time-shift each one of the curves such that their initial 
values would sit on the 60°C — > 30°C curve. The failure 
of the time-shifted curves to collapse on a master curve 
implies that using a single non-equilibrium variable (here 
the deviation of the volume from its equilibrium value) 
would be inadequate for constructing a predictive model 
of structural relaxation. 
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FIG. 1. (color online) (a) Down-jumps with a fixed target 
temperature T = 30°C in PVAc. Data digitized from Fig. 16 
in [lj]. (b) The above measurements time-shifted onto the 
60°C -> 30° C curve. 



This observation seems to be in line with the com- 
mon knowledge that glassy relaxation is characterized 
by a broad spectrum of relaxation times. A fundamen- 
tal modeling approach would incorporate these various 
timcscales in the time evolution of the probability distri- 
bution function of the volume v of mesoscopic material 
elements f(v,t), as it approaches the stationary equilib- 
rium distribution f^ q (v) at T during structural relax- 
ation. This is a dauntin g ta sk that has been pursued 
only in simple models [20I. J2ll]. Our goal here is to show 
that while models that use only the macroscopic volume 
V(t) = N J^°vf(v,t)dv (where N is the number of ele- 
ments) cannot be complete, they can still teach us some- 
thing and might serve as a starting point for constructing 
an adequate phenomenological model. 



There are two basic approaches for understanding log- 
arithmic relaxations. The first approach sees the relax- 
ation as a linear response, i.e. relaxation rates are in- 
dependent of the state of the system. A logarithmic re- 
sponse can then be obtained by summing over a spec- 
trum of exponential relaxation modes. This approach 
was suggested in jlj], pursued in 22 1, and was recently 
also invoked in the context of electron glasses Es- 
sentially, the evolution of the deviation from equilibrium 
is assumed to take the form 8(t) oc J °°p(A) exp(— Xt)dX, 
where p(X) is the distribution of relaxation rates. When 
p(X) ~ A^ 1 in a certain range, logarithmic behavior 
emerges. The second approach, suggested in many con- 
texts (e.g. in Q), describes logarithmic relaxation as a 
result of the dependence of rate on instantaneous state, 
e.g. slowing-down due to compaction. A nonlinear equa- 
tion of the form <5oc — cxp(65) is then proposed to yield 
£oc-log(i)/6. 

We follow the latter approach and propose to describe 
structural relaxation by the following equation 



5 = a (1 — e 



b 6\ 



(1) 



Here a is a basic relaxation rate and & is a constant. 
The 1 added in brackets ensures that 8 vanishes with 8. 
Equation ([!]) admits the analytic solution 
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S(t) = --log [1-(1- 



(2) 



where <5o is the initial condition. For large enough initial 
amplitudes, for which exp(— bSo) <C 1, and short times 
abt <C 1, we have 8(t) ~ — \og{abt)/b. The final stage of 
the relaxation is exponential, S(t) ocexp(— abt). For large 
negative initial amplitudes, for which exp(— b8o) 3> 1, 
and short times at <C — <5oi the relaxation is linear 8 ~ 
So+a t. Again, the final stage of relaxation is exponential, 
8(t) oc — exp(— abt). This marked asymmetry between 
8 > and 8 < relaxations naturally emerges from the 
exponential 8— ¥— 8 asymmetry in Eq. ([lj. 

To test the model we present in Fig. [5] fits of Eq. ^ 
to a significant portion of Kovacs' original data (l4l| . For 
each curve the parameters 8q, a, b were independently var- 
ied. Note, however, that <5o is essentially determined by 
the first data point and hence is not a real fitting parame- 
ter. Down-jumps in PVAc with fixed target temperature 
T (panel (a)) and fixed initial condition To (panel (b)) 
are satisfactorily captured by Eq. We checked the 

equation also against data for glucose (panel (c)) which 
appears in [lij]. Panels (d-e) show fits to up-jumps - ex- 
periments where the initial temperature To is lower than 
T, making the initial condition 8q < 0. Panel (d) is 
plotted in a linear scale to show the manifestly linear 
portion of the relaxation. Usually both up-jumps and 
down-jumps are plotted using a logarithmic time scale to 
show the asymmetry between the two responses. Here 
we suggest that this asymmetry can be understood to 
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FIG. 2. (color online) Fits of Eq. © to Kovacs' data [Tj. The solid black curves are three-parameter fits, (a) Down-jumps in 
PVAc with fixed T = 30°C. See Fig. □ (b) Down-jumps in PVAc from the same initial condition T o = 40°C to T = 25, 37.5° C 
(top to bottom), (c) Down-jumps in glucose from a fixed initial temperature To = 40°C to T= 19.8, 32.5°C (top to bottom), 
(d) Up-jumps in PVAc. Note the linear scale. Here Tb = 30, .., 37.5°C (bottom to top) to T — 40°C. The measurement starts at 
ti — 36 sec. (e) Up-jumps in glucose where Tb = 25°C is fixed while T = 30, ..,37.4°C (bottom to top). The measurement starts 
at U = 36 sec. 



emerge from the nonlincarity of the response. Additional 
data were shown to be consistent with Eq. ([I]) 23 1. 

Is it a mere curve fitting? Equation {1} provides an ex- 
cellent description of the data, at a price of independently 
varying the model parameters for each set of experimen- 
tal conditions. In addition, we already know that a single 
variable approach as in Eq. ([lj cannot constitute a pre- 
dictive model. Nevertheless, we argue that the variation 
of the parameters a and b with experimental condition 
may be physically meaningful and provide us with phys- 
ical insight into structural relaxation. Indeed, Fig. [3] 
shows that a and b do show systematic variations with 
To and T. Panel (a) shows a(T) for PVAc and glucose 
under a fixed initial temperature To (in each experiment, 
both down and up-jumps). The first observation is that 
a(T) is a function, i.e. it exhibits a smooth variation 
with T. The same is true for b(T), see panel (b) for 
measurements with fixed Tq. In fact, a exhibits a strong 
exponential dependence on T. This is a manifestation of 
the dramatic slowing down of dynamics associated with 
the glass transition. The possible origin of this depen- 
dence in thermally activated processes will be discussed 
below. 

In the present context, the observation that the struc- 
tural state of a glass cannot be described by a single non- 
equilibrium variable implies that a and b should depend 
on the initial temperature Tq. Indeed panel (b) demon- 



strates that b depends on To. Here the parameter b in 
PVAc and glucose is plotted against the jump size T — T 
for both down-jumps and up-jumps and a wide range of 
variation of both T and Tq. This shows that b is in fact a 
function of both T and To, which strikingly implies that 
a relaxing glass carries memory of its original state at To 
for very long times, sometimes for days! The dependence 
of a on initial To (not shown here) is the opposite of that 
of b. In spite of this fact, we have not found any sim- 
ple connection between the two. We also verified (not 
shown) that the final exponential relaxation rate, which 
is controlled by the product a(T, To) 6(T, To), depends on 
T) as well. 

Another interesting feature of panel (b) is the approx- 
imate collapse of all measurements on a single curve as 
a function of To — T. This lends support to the view 
that b has a physical meaning. We also note the ap- 
parent discontinuity of b when passing from up-jumps 
to down-jumps. Discontinuity of this kind in relaxation 
times is known in the literature as the T e f /-paradox (3] ■ 
We take it to mean that as the asymptotic volume at 
T is approached from above (through a down-jump) or 
from below (through an up-jump), the system explores 
disparate regions in phase space. This is yet another in- 
dication that the volume alone does not tell the whole 
story about glassy relaxation. This might also be related 
to the fact that the volume and other thcrmo-mechanical 
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FIG. 3. (color online) Variation of the parameters in Eq. 
(JXJ) with initial temperature To and target temperature T. 
Legend: (□) PVAc down-jumps from T = 40° C [3 (x) 
PVAc down-jumps to T = 30°C (A) PVAc Thumps to 
T — 40° C H (+) PVAc T-jumps to T — 35° C [jj, H| (o) 
PVAc down-jumps from T = 35°C [H (0) PVAc down- 
jumps from To = 40° C [28|] (□) Glucose down-jumps from 
T = 40°C* (A) Glucose up-jumps from T = 25°C [3 

(a) Variation of relaxation rate a with T for experiments of 

fixed initial condition. The dashed line is an exponent e T ' T 
with T — 1°K. The error bars are smaller than the symbols. 

(b) For PVAc and glucose b shows a consistent dependence 
on jump magnitude To — T. 



properties of glasses do not necessarily equilibrate simul- 
taneously (see [l[ section B.1.8). 

Equation (JT|), with its parameters a and 6, does seem to 
contain meaningful physical information about the pro- 
cess of structural relaxation. In order to quantify this 
information, we propose to rationalize the equation us- 
ing a simple two-state model. A common way to moti- 
vate such an equation is to assume that an Arrhcnius-likc 
process depends on an observable. In this spirit, let us 
assume that the volume evolves through activated jumps 
of volume elements between a contracted "— state" and 
an expanded "+ state" . A rate equation that describes 
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FIG. 4. (color online) The volume dependent energy barrier 
E-(Vt) = — fcsT log(a/fco) for the parameter a in Fig. [3] 
panel (a). Here kg = 10 12 sec -1 . Solid lines are linear fits. 



this evolution is 



h = —k—Ti + fc+(l — n) , 



(3) 



where n is the fraction of elements in the + state. A 
key assumption here is that the activation energies are 
volume dependent, i.e. k± cx exp(— E±(n)/kBT). E±(n) 
are energy barriers for expansion and contraction. They 
are assumed to grow with density - it is harder to move in 
a denser system. Taking into account the fact that 6V/V 
is small, of the order of 10 -3 , we approximate n~riT+6n, 
where ut is the equilibrium value of n at temperature T. 
By neglecting terms of order Sri, using the equilibrium 
condition k-(jiT)riT = fc+(nx)(l — n?) and writing 5 n = 
Sn/riT, we find that 



S„ = k-{n T ) (e 



A+<5„ 



„A-6„ 



(4) 



where A± = — (ht / 'kBT)dE±(riT) / 'dn are expected to be 
non-negative, as barriers should become smaller with ex- 
pansion. Equation ^ then emerges as a special case 
when a(T) = k-(n T ), A + = and b(T)=A_. 

This model allows us to extract an energy barrier scale 
from a = k n cxp(—E-(Vr)/kBT), where Vr is the equi- 
librium volume at T. Since fco is a pre-exponential fac- 
tor, its exact value is not important. We take fco ~ 10 12 
sec -1 , which is the typical scale of molecular vibration. 
In Fig. [4] we extract values of the energy barrier E_ 
from the experiments summarized in Fig. OHa). The 
activation barrier turns out to be on a leV scale, which 
might be indicative of cooperative rearrangements of tens 
of monomers. 

To further test the latter possibility, we estimate the 
size of rearranging regions by a dimensional argument 
using the measured derivative dE-(Vr)/dT. We define 
a volume scale v as 



nV 



dE_ 



ndE- 



a dT 



(5) 
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where k = V~ 1 dV/dP is the isothermal compressibility 
and a = V~ 1 dV/dT is the thermal expansion coefficient. 
Linear regressions shown in Fig. @]yield dE_ / dT ~ —0.02 
cV/K. Taking for PVAc a ~ 7x 10" 4 K" 1 and k~0.42 
GPa -1 [3l we find v ~ 2 nm 3 . A similar estimate for 
glucose also results in a cubic nanometer scale. These 
estimates imply that structural rearrangements involve 
tens of basic units (e.g. monomers in a polymer), which 
seems consistent with the typical size of dynamical het- 
erogeneities reported in the literature [29(. Therefore, 
while the model cannot explain the dependence of the 
parameters on the initial temperature To, it does seem 
to be sensitive to the dominant scales of the underlying 
relaxation processes. 

In summary, we proposed a simple mean-field equation 
to describe structural relaxation in glasses. Our analysis 
clearly demonstrates that the structural state of a re- 
laxing glass cannot be fully accounted for using a single 
non-equilibrium variable. In fact, we explicitly showed 
that a glass may carry information about its thermal 
history for extremely long timescales. Furthermore, the 
simple analysis allows the estimation of typical energy 
and volume scales associated with thermally-activated 
structural relaxation processes. Theses estimates indi- 
cate the existence of cooperative rearrangements on a 
super-molecular scale. 
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